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Abstract  of  Dissertation  Presented  to  the  Graduate  School 
of  the  University  of  Florida  in  Partial  Fulfillment  of  the 
Requirements  for  the  Degree  of  Doctor  of  Philosophy 

GRAND  UNIFIED  SUPERSYMMETRIC 
INFLATIONARY  COSMOLOGY 

By 

Marti  Ruiz-Altaba 
May  1987 

Chairman:  Pierre  Ramond 
Major  Department:  Physics 

A  theoretical  laboratory  is  presented,  within  the  framework  of  supergrav- 
ity  unification,  which  satisfies  all  experimental  and  cosmological  constraints 
below  the  Planck  scale,  notably  those  from  standard  big  bang  nucleosynthe- 
sis, baryogenesis,  and  the  fermionic  spectrum.  An  intermediate  scale,  around 
1010  -  10nGeV,  arises  from  the  examination  of  inflation,  supersymmetry 
breaking,  the  Peccei-Quinn  symmetry,  fermion  masses,  proton  decay,  and  the 
baryon  asymmetry  of  the  universe.  Careful  attention  is  paid  to  the  renormal- 
ization  of  proton  decay  operators.  Although  a  realistic  low  energy  spectrum 
follows  from  an  isospin-like  global  symmetry  (familiarity),  a  rather  baroque 
higgs  sector  is  needed  to  cancel  anomalies  and  implement  the  missing  mul- 
tiplet  mechanism.  A  low  reheat  temperature  and  baryogenesis  from  out-of- 
equilibrium  decays  are  perhaps  the  most  noteworthy  aspects  of  the  model, 
whose  general  success  is  not  yet  matched  by  any  models  derived  from  strings 
rather  than  point  fields. 


v 


INTRODUCTION 

One  of  the  most  astounding  successes  of  the  grand  unification  program,1 
initiated  about  ten  years  ago,  was  the  understanding  of  the  origin  of  the  baryon 
asymmetry  of  the  universe.2  The  program  suffered  an  early  setback  when  min- 
imal SU(5)  failed  to  explain  quantitatively  this  baryon  asymmetry.  It  also 
predicted  a  fast  proton  decay  rate  which  was  soon  ruled  out  by  experiment. 
Supersymmetry3,4,5  generally  extends  the  proton's  lifetime  and  prefers  the 
strange  and  muonic  decay  modes.  The  experimental  limits  on  the  proton's 
lifetime6  and  the  requirements  of  efficient  baryogenesis  tightly  constrain  the 
parameters  in  the  Lagrangian  of  candidate  unification  models.  An  excellent 
theoretical  laboratory  arises  from  the  locally  supersymmetric  grand  unified 
model  below,7'8'9,10'11, 12,13  which  is  consistent 

a)  internally  below  the  Planck  scale  (anomaly-free), 

b)  with  the  standard  SU(3)  x  SU(2)  x  U(l)  model  at  accelerator  energies, 

c)  with  big  bang  inflationary  cosmology14,15'16  at  very  early  times,  and 

d)  with  underground  detector  limits  on  proton  decay. 

The  general  traits  of  this  supersymmetric  inflationary  cosmology  (SIC)  are 
presented  below.  The  seminal  papers  for  SIC  are  Refs.  7  and  8.  Further  work 
appeared  in  Refs.  11  and  12.  Although  the  model  is  not  particularly  aesthetic, 
it  serves  to  prove,  by  construction,  the  consistency  of  all  known  constraints 
with  simple  supergravity  unification. 

Supersymmetry   is   spontaneously    broken    by    the  O'Raifeartaigh 
mechanism17'18  and  its  breaking  is  related  to  the  end  of  the  inflationary  period 
in  the  evolution  of  our  universe.  The  cosmological  constant  can  be  set  to  zero 


while  supersymmetry  is  broken,  but  a  contribution  of  0{M^Mpl&nciL)  cannot 
be  avoided  after  the  observable  sector  finds  out  about  supersymmetry  break- 
ing. No  mechanism  is  known  in  any  model  to  solve  the  cosmological  constant 
problem,  so  we  proceed  with  the  rest  of  the  program  without  much  uneasiness. 

A  specific  model  with  SU(5)  invariance  is  written  down,  relying  freely  on 
insights  gathered  from  explorations  with  larger  groups.  The  mass  spectrum  of 
leptons  and  quarks  is  incorporated  at  the  tree  level  with  the  aid  of  couplings 
between  matter  and  a  nonminimal  higgs  sector  including19  a  45.  The  element 
KMcb  of  the  quark  mixing  matrix  has  been  measured  from  the  (long)  B- 
meson  lifetime.20  We  incorporate  it  and  a  low  top-quark  mass  (claimed  and 
subsequently  disclaimed  by  the  UA1  group21),  obtaining  reasonable  fermion 
mass  relationships  and  adequate  CP  violation. 

Since  gravitinos  are  created  during  reheat,22  the  temperature  Tr  of  the 
heat  bath  generating  them  is  bounded  above.23  If  reheated  fields  are  thermal- 
ized  immediately,  the  products  of  the  reheating  process  are  restricted  to  having 
masses  <  Tr,  lest  their  abundance  be  strongly  suppressed,  ~  exp(— Tji/m). 
However,  if  the  reheat  is  not  an  equilibrium  process,24  then  no  such  strict 
bound  can  be  derived  for  the  reheated  products,  which  could  be  more  massive 
than  Tr.  In  the  SIC  theory  presented  here,  reheat  is  caused  by  the  decay  of 
inflatons,  of  mass  around  the  intermediate  scale,  1010  —  1011  GeV,  and  yet  the 
reheat  temperature  is  only  about  106  GeV. 

Higgs  triplets  of  masses  around  the  intermediate  scale  can  be  used  for 
generating  the  baryon-number  asymmetry25'26'27  without  causing  fast  proton 
decay  via  the  dimension  five  higgsino-mediated  operators.28'29'30  The  method 
presented  here  calls  for  right-handed  neutrinos  with  masses  around,  again,  the 
intermediate  scale. 


Proton  decay  rates,  with  largest  branching  ratios  into  strange  and  muonic 
modes,  put  a  lower  bound  on  the  masses  of  the  higgs  triplets  whose  interactions 
violate  baryon  number.  The  requirement  that  these  same  triplets  be  produced 
from  inflaton  decay  to  generate  the  baryon  asymmetry  of  the  universe  sets  an 
upper  bound  on  their  masses.  A  narrow  window  satisfying  these  two  bounds 
is  found  at  the  intermediate  scale. 

This  dissertation  is  organized  as  follows.  First,  the  SIC  superpotential  is 
presented  and  the  hidden  sectors  are  glossed  over.  The  grand-unified  (GUT) 
sector  is  then  presented  in  detail  along  with  its  implications  for  fermion  masses, 
CP  violation,  and  neutrino  oscillations,31  in  terms  of  the  Yukawa  couplings  of 
the  superpotential  and  the  various  symmetry-breaking  scales.  We  then  com- 
pute the  decay  rates  for  the  proton,  under  flavor-helicity  SU (6),  paying  careful 
attention  to  the  renormalization  of  the  effective  B-violating  low-energy  opera- 
tors. Finally,  we  study  baryogenesis  and  present  a  mechanism  involving  heavy 
right-handed  neutrinos,  whose  masses  are  also  around  the  intermediate  scale. 

The  field  content  of  supergravity32,33  is  as  follows.  Chiral  superfields 
(Da<f>{  =  0)  contain  matter,  in  spinorial  and  scalar  varieties.  Gauge  super- 
fields  (Wa  —  D  DaV,  V  =  V>)  contain  gauge  vector  bosons  and  the  spinorial 
gauginos.  In  addition,  the  spectrum  contains  the  gravity  supermultiplet,  with 
the  graviton,  the  gravitino  and  the  Kalb-Ramond  field.  Supergravity  theories 
are  never  torsion-free  due  to  the  interdependence  of  the  fermionic  and  bosonic 
coordinates  of  superspace.34 

Minimal  N  =  1  supergravity  in  four  dimensions  has  two  arbitrary  indepen- 
dent functions,  which  reduce  to  only  one  superpotential  P  in  the  "minimal" 


case.  Kinetic  terms  for  all  fields  are  in  this  case  canonical,  and  the  only  re- 
quirement on  P  is  that  it  be  analytic,  i.e.,  function  only  of  the  superfields  and 
not  of  their  conjugate  partners. 

The  superpotential  for  the  SIC  model  consists  of  three  sectors, 

P  =  I  +  G  +  S  (1.1) 

where  fields  in  one  sector  do  not  appear  in  the  other  two.  Both  I  and  S  contain 
only  gauge  singlets,  responsible  for  inflationary  expansion  and  supersymmetry 
breaking  respectively.  The  sector  G  contains  all  the  usual  matter,  the  higgs 
sector,  and  more. 

The  purely  scalar  potential  for  a  minimal  supergravity  theory  with  super- 
potential  P  is  given  by 
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scalar  exP 


— ipi2 


+  D  —  terms 
(1.2) 


where  the  sums  run  over  all  scalar  (super) fields  in  the  theory,  denoted  symbol- 
ically by  fa,  and  M  is  the  reduced  Planck  mass 

M  =  -±=MPlanck  =  2.4xl018GeV  (1.3) 

V87T 

Cross-talk  between  the  sectors  is  always  suppressed  by  powers  of  M.  Each 
of  the  three  sectors  interacts  with  the  other  two  only  gravitationally,  so  thermal 
equilibrium  is  possible  only  within  each  sector  or  as  a  consequence  of  fine-tuned 
initial  conditions. 


INFLATION 


The  inflation  sector  J  of  the  superpotential  can  be  constructed  very  simply7 
by  assuming  it  contains  only  one  chiral  superfield,  $.  Scaling  all  fields  by  M 
and  considering  the  scalar  potential  due  to  /  alone, 


Vj 


1*1' 


-  31/1 


(2.1) 


At  the  minimum  of  this  potential  (the  true  vacuum  up  to  interference  with 
other  sectors),  the  cosmological  constant  should  vanish,  so  Vj($q)  =  0.  Re- 
quiring further  supersymmetry  to  be  unbroken  at  this  minimum  implies  that 
the  gravitino  mass  must  be  zero,  which  imples  that 


91  *tr 
^  +  $  1 


=  0 


$=$o 


Hence 


mo)  =  ||(*o)  =  0 


(2.2) 


(2.3) 


A  very  simple  J  satisfying  the  above  requirements  is  a  quadratic  superpo- 
tential 

J  =  A2($-$0)2  (2.4) 

where  A  is  a  mass  parameter  setting  the  scale  for  inflation.  This  is  the  small 
parameter  which  appears  in  all  succesful  inflationary  models,2  fixed  around 
10-4M  to  avoid  the  supersymmetry  entropy  crisis  or  Polonyi  problem.8'22 
Without  loss  of  generality,  we  can  take  the  vacuum  value  of  the  inflaton, 
(<£)  =  $o>  to  be  real. 
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Requiring  the  scalar  potential  to  be  flat  at  the  origin, 


( 


d<f>. 


(0)  =  0  (2.5) 


so  that  (new)  inflation  can  occur  with  a  slow  rollover,  sets  $o  =  M  =  1. 
Letting  the  scalar  component  of  the  chiral  superfield  $  be  <f>e%6 the  scalar 
potential  from  this  sector  only  is 

Vj  =  A  V2  <j>6  -  4(f>5  cos  0  +  (3  +  4  cos2  0)<t>4  -  4  cos  0<f>3  +  (3  -  4  cos2  0)<f>2  +  1 

(2.6) 

For  any  value  of  <p  ("the  radius")  this  potential  has  a  minimum  at  0  =  0.  It  is 
the  only  local  minimum,  so  there  is  no  ^-degeneracy  and  no  domain  walls  may 


arise. 


W)  =  W,  8  =  0)  =  A  V2  (<£6  -  4<f>5  +  7<j>4  -  4<j>3  -  <j>2  +  l) 

=  A4  (l  -  4*3  +  j <t>4  -  *<f>5  +  | **  +  <W7)) 

=  4eA\<f>  -  l)2  +  12eA4(<£  -  l)3  +  0{{<f>  -  l)4) 

(2.7) 

whence  the  mass  of  the  inflaton  is 

=  2v/eA2  (2.8) 

Similarly, 

Vj{0)  =  VI{0,<t>  =  l)  =  8eA4(l-cos0)  (2.9) 
so  the  mass  of  the  angular  field  0  is 

me   -    y/im^  =  2\/2eA2  (2.10) 

The  potential  (2.7)  is  a  standard  inflationary  potential,  flat  at  the  origin, 
with  a  slow  descent  towards  the  minimum  at  M,  and  a  steep  divergence  for 
values  of  the  field  larger  than  M .   The  inflaton  <f>  has  cubic  self-interaction 
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of  strength  12A2,  and  linear  (gravitational)  interactions  with  other  sectors  of 
strength  ~  A2  as  well.  The  inflaton  will  couple  to  fields  in  S  and  G,  with  a 
preference  for  the  heaviest  ones,  and  a  decay  rate 

I>  «  (A2)2m^  cA6  (2.11) 

At  Planck  time  all  particle  species  in  the  universe  are  maintained  in  thermal 
equilibrium  by  gravity.  Regardless,  however,  of  the  initial  conditions,  each  of 
the  three  sectors  decouples  from  the  others  while  gravity  becomes  weaker  as 
the  universe  expands.  Fields  within  each  sector  thermalize  through  Yukawa 
and  gauge  interactions,  so  the  inflaton  superfield  effectively  evolves  alone  in  a 
thermal  bath  provided  by  S  and  G. 

The  inflaton  will  roll  away35  from  its  initial  configuration36'37 

<jxH0,    <j>  ~  0  (2.12) 

destabilized  by  the  (otherwise  negligible)  Hawking  radiation  term 

AVH  ~  -AV  (2.13) 

As  the  universe  expands,  the  energy  density  becomes  dominated  by  the  po- 
tential energy  of  the  inflaton,  which  can  reasonably  be  assumed38  to  evolve 
according  to  the  classical  equations  of  motion 

j>  +  3H<j>+  =  0  (2.14) 

a<p 

"2  =  ip  (I*2  +  v'<*>)  <215> 

where  H  is  the  Hubble  constant  (parameter,  rather)  with  its  value  in  anti- 
deSitter  spacetime  given  by 
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Inflation  occurs  with  $  «  3H0  <$>  from  a  time  t*  when  <f>  ~  H0  until  a  time 
tj  ~  mT1.  The  inflaton  then  starts  evolving  like  matter,  until  a  time  i#  ~  T^1 
when  it  decays  into  radiation  and  reheats  the  universe.7,8  Using  conservation 
of  energy,  the  increase  in  radiation  energy  density  is  determined  by 


P<f>{tR)  =  P+{tj) 


2 


=  A4 

~  A12 

TT2 

*TR 

(2.17) 


with  g*  the  effective  number  of  light  degrees  of  freedom  at  the  time  of  reheat 
tft.  The  reheat  temperature  Tr  for  the  heat  bath  in  the  causal  bubble  is,  by 
definition, 

TR  «  A3  (2.18) 

To  avoid  the  monopole  problem,14,39  must  be  less  than  the  unification 
scale=  Mx\  hence  A2  <  10~2.  Below  we  shall  find  that  actually  A  ~  10-4, 
so  the  reheat  temperature  Tr  is  only  about  106GeV.  Despite  this  low  reheat 
temperature  (which  does  not,  by  far,  produce  any  problems  with  nucleosynthe- 
sis), inflatons  are  produced  copiously  and  inflaton  decay  does  generate  enough 
baryon  asymmetry,  due  to  the  lack  of  equilibrium  among  the  various  hidden 
sectors. 

The  number  Nj  of  e-folds  of  inflationary  expansion  (Rr  ~  eNlRj)  will  be9 

f<f>reheat 

Nj  =    /  H0dt  (2.19) 

JfcnitiaX 

This  is  easily  much  larger  than  the  observationally  required  Nj  >  65.  Nj  is 
around  108  if  fa  <  Hq  ,  <f>r  ~  M/10,  so  there  is  a  large  set  of  initial  conditions 
leading  to  sufficient  inflation. 
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In  the  deSitter  (inflationary)  stage,  the  quantum  fluctuations40,41,42,43  of 
the  inflaton  are  of  order  Hq.  Their  comoving  size  leaves  the  horizon  during  in- 
flation (Hubble  radius  is  constant),  and  they  reenter  the  horizon  after  inflation, 
thus  generating  a  spectrum  of  density  fluctuations  at  scale  A: 

2 

(2.20) 

Approximating  Vj  by  its  truncated  Taylor  series  around  the  origin,  Eqs. 
(2.14)  and  (2.15)  can  be  solved  for  <j>{tj)  to  give7 

^  -  104A2  ~  1(T4  (2.21) 
P 

which  is  the  preferred  value  for  galaxy  formation  with  adiabatic  density  fluc- 
tuations, known  in  the  trade  as  the  "cold  dark  matter  scenario."  Interestingly, 
Hawking  has  shown44  that  if  <  10-4,  then  the  mass  of  the  inflaton  must 
be  less  than  1015  GeV,  or  A  <  10~1,5.  Let  us  emphasize  that,  given  the  value 
of  A  ~  10-4  required  from  the  discussion  of  the  gravitino  problem,8  the  mag- 
nitude of  the  density  fluctuations  comes  out  as  a  prediction  of  the  theory. 


4>Htr) 


1  +  -log(Atfo) 
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SUPERSYMMETRY  BREAKING 

Consider  now  the  supersymmetry-breaking  piece  S  of  the  superpotential 
P.  It  cannot  be22  of  the  Polonyi45  type  (linear  in  the  superfield).  Take  S  to 
be  of  the  O'Raifeartaigh  type,8'17 

S  =  A{\B2  +  fi2)  +  mBC  +  n  (3.1) 

where  A,  B,  C  are  chiral  superfields  and  A,  /x,  m,  and  n  are  parameters.  We 
can  set  m  ~  \i  and  adjust  n  to  ensure  that  the  cosmological  constant  vanishes 
at  the  minimum  of  the  full  potential,  taking  into  account  the  inflation  and 
gauge  sectors. 

The  F  term  §4  breaks  supersymmetry17'18  so  that  the  spinor  compo- 
nent of  A  is  the  Goldstone  fermion  which  provides  the  gravitino  with  the 
helicity  ^  states  it  needs  to  become  massive.  The  scalar  component  of  A, 
the  O'Raifearton,  couples  weakly  and  stores  potential  energy  during  inflation. 
Then  it  evolves  like  matter  until  it  decays  preferentially  into  a  gravitino  pair 
through  a  one-loop  diagram  with  three  internal  B  propagators,  with  rate 

1^3  3  m  TtVm  (3.2) 

where      =  A2/47r,  and  a\n  is  the  one-loop  O'Raifearton  mass. 

The  gravitino  mass  is  mz  =  H2/M,  and  the  B,  C  fields  have  a  mass  of  order 

3 

fj,.  For  supersymmetry  to  solve  the  gauge  hierarchy  problem46'47,48,49'50,51  (its 
precise  raxson  d'etre),  mz  must  be  at  most  around  a  TeV  so  /x  ~  1010-11  GeV 

2 

~  (10-7  -  10_8)M. 
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Now  gravitinos  couple  to  matter  very  weakly  (with  gravitational  strength), 
so  they  are  long-lived  particles.  They  may  thus  decay  into  photons  and  photi- 
nos  (lightest  odd-fl  particles,  R  =  (_i)3^+^+25)  after  the  onset  of  nucle- 
osynthesis. Their  number  density  at  decay  is  severely  constrained52,53,54  by 
the  success  of  conventional  nucleosynthesis55  in  explaining  the  primeval  abun- 
dances of  light  element.  The  decay  of  gravitinos  into  gluon-gluino  pairs,  for 
instance,  can  dramatically  alter  the  deuterium  abundance.56,57 

The  abundance  of  gravitinos  at  the  time  of  nucleosynthesis  is  constrained, 
whence  the  reheating  temperature  is  bounded24 

1  TeV' 


TR  <  2.5  x  107  GeV 
Using  Tr  ~  A3  and  mz  ~  /z2,  this  implies 


f/l3 

2 


(3.3) 


A  V  <  4.2  x  10"27 


(3.4) 


which  entails,  under  the  requirement  (for  the  solution  of  the  hierarchy  problem) 
that  n  ~  10~8,  that8 


A  <  1(T4  (3.5) 
If  gravitinos  happened  to  be  stable,  then  the  bound  (3.3)  is  less  restrictive,24 


TR  <  1.3  x  1011  GeV 


1  TeV 

T7l3 

2 


implying 


A  <  10 


-2.5 


(3.6) 


(3.7) 


To  make  sure  that  infiaton  production  of  gravitinos  is  small,  since  the  infla- 
ton  couples  most  strongly  (~  A2)  to  the  heaviest  fields,  we  must  kinematically 
forbid  the  infiaton  from  decaying  into  O'Raifeartons,  as  well  as  B  and  C  fields, 


m<f>  <  mA->  mB  >mC 


(3.8) 
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and  require  also  to  be  small,  <  10-3,  in  order  to  suppress  one-loop 
effects.8 

A  careful  analysis8  of  these  one-loop  effects  reveals  a  lower  bound  for  the 
inflation  scale 

A  <  10"4-2  (3.9) 
whence  A  is  rather  tightly  constrained,  much  to  nature's  credit. 
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GRAND  UNIFIED  SECTOR 


11 


(4.1) 


We  write  the  superpotential  for  matter  fields  as 

„  /  0     a'H     0  \  „  /  0      dK       0  \ 

Gy  =  Fr     oH    cL     0   J  T  +  T     I  dK    eK'    fK"  T 
V  0       0     bHj  V  0     fK"     gK  J 

(   0      dxK       0  \ 
+  ST     d2K     eK'     hK"    F  +  {j^S^  +  y25353) 

V  o  /2x"  sir  7 

where  TT  =  (ri,r2,r3),  F  =  {Fl,F2,F3),  and  Sr  =  (SUS2S3)  are  the  N  = 
1  chiral  superfields  containing  matter  in  three  generations  or  families,  and 
H,  L,  K,  K'  and  K"  are  higgs  superfields,  and  4>\  is  a  group  singlet  whose 
vacuum  expectation  value  (<£i)  =  Mq  =  1010GeV  breaks  the  Peccei-Quinn 
symmetry58. 

The  explicit  field  content  of  Fn,  Tn,  and  Sn  is  shown  below,  where  Greek 
indices  run  from  one  to  three  and  denote  color.  All  fields  are  chiral  super- 
fields  with  left-handed  Weyl  spinor  components.  The  fields  Fn  are  vectors,  Tn 
antisymmetric  tensors,  and  Sn  singlets: 

Fn  =  tn 


1     (ZafaiO1     -«na  ~4 
Tn  =  ~E  [  0 


v/2 


na 

c 


0 


(4.2) 


Sn  =Vn 

The  representations  and  global  symmetry  charges  of  the  various  fields  in 
G\  are  shown  in  Table  I.  The  only  unusual  feature  is  the  higgs  field  L,  which 


14 

transforms  like  a  45.  The  first  two  terms  in  (4.1)  constitute  the  Georgi-Jarlskog 
model19  with  an  extra  eK'T2T2  coupling,  needed  to  correct  the  prediction 
KMci  =52  incompatible  with  the  value  inferred  from  measurements  of  the 

B-meson  lifetime20  [KMcb  =  0.0435  ±  0.0003)  and  the  small  top-quark  mass 
(<  150  GeV)  required  by  unitarity. 


TABLE  I.  Field  content  and  symmetries 
orial,  Peccei-Quinn,  familiarity)  of  G\. 

(vect- 

Superfield 

SU(5)  irrep 

V 

X 

7 

5 

-3 

1 

ri 

10 

1 

1 

Si 

1 

5 

1 

F2 

5 

-3 

1 

T2 

10 

1 

1 

S2 

1 

5 

1 

1  3 

^ 

o 

xJ 

0 

r3 

10 

1 

0 

S3 

1 

5 

0 

H 

5 

2 

-2 

0 

L 

45 

2 

-2 

2 

K 

5 

-2 

-2 

0 

K' 

5 

-2 

-2 

2 

K" 

5 

-2 

-2 

1 

4>1 

1 

-10 

-2 

0 

The  four  phase  symmetries  can  be  understood  as  follows.11  First,  there  is 
the  vectorial  symmetry  V,  which  becomes  B—L  below  the  electroweak  breaking 
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scale.  It  forbids  unwanted  terms  like  FK  and  FFT  and  has  no  anomalies. 
Second,  we  have  the  anomalous  Peccei-Quinn  symmetry  X,  needed  to  solve  the 
strong  CP  problem.59  Third  is  the  chiral  R  symmetry  of  cubic  superpotentials, 
broken  by  gravitational  couplings  to  the  hidden  sectors  /  and  S.  Lastly,  we 
have  the  J  symmetry,  familiarity,  which  is  the  only  one  to  distinguish  among 
the  three  families  or  generations.  Familiarity  7  is  like  family  isospin,  traceless 
in  family  space. 

Crucially,  we  forbid11  mixing  terms  between  the  two  sets  of  higgs  fields 
with  different  V,  such  as  HK,  HK',  HK"  or  LEff,  LEif',  LEK",  where 
E  ~  24  or  75.  Such  mass  terms  would  lead  to  disastrously  fast  proton  decay 
rates,  through  dimension-five  operators  (see  the  discussion  below  on  proton 
decay) . 

As  a  further  and  rather  innocent  simplification,  we  shall  assume  that  the 
matrices  are  symmetric,  i.e.,  that  o  -  a',  d\  —  d>2,  and  f\  =  fa. 

Fermion  Masses 

The  higgs  superfields  produce  after  SU(5)  breaking  light  (~  M\y)  doublets 
and  heavy  (~  M^/M)  triplets,  via  the  missing  multiplet  mechanism.60'61  Our 
purpose  now  is  to  determine  the  fermion  masses,  in  order  to  constrain  the 
various  Yukawa  parameters  appearing  in  the  superpotential.  The  A/jy  =  | 
masses  come  from  the  5^(2)^  x  U(l)y  — >  ^(l)em  phase  transition,  when 
neutral  components  of  higgs  SU{2)i  doublets  acquire  a  vacuum  expectation 
value  (vev)  a.  The  Alyy  =  1  masses  come  from  the  breaking  of  the  Peccei- 
Quinn  phase  symmetry,  when  some  singlets  acquire  vet/'s  Mq.  Rewriting 
(4.1)  in  terms  of  the  higgs  (3,1)  triplets  and  (1,2)  doublets,  we  arrive  at  the 
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SU(3)  x  SU(2)i  x  U(l)y  superpotential  for  this  sector,  with  due  care  for  the 
normalization  of  the  kinetic  terms  of  the  component  fields. 

The  notation  ipx  is  shorthand  for  tpTCx,  with  C  the  charge-conjugation 
matrix  and  T  the  transposition  in  Dirac  space.  A  useful  identity  is 

($L)    C  =  rPR  (4.3) 

Since  all  spinors  are  left-handed, 

rf>TCX  =  rPlCxL  =  *RXL  =  XlC1>L  =  X°R^L  (4-4) 
In  the  more  convenient  two-component  notation,  in  which 

Vc  -  o2tfR  (4.5) 
the  various  products  are  of  one  of  the  forms  below: 

V>V  -  M)T(*2XR)  =  -<l>RxR  =  -Wrxr)* 

(4,6) 

V»CX  ~*  {^*r)  XL  =  -1>R?2XL  =  -$RXL 

<1>XC  ~*  ^l(°2Xr)  =  -{1>\°2XRY  =  -$LXR)* 
corresponding  to  the  left  Majorana  mass  term,  the  right  Majorana  mass  term, 
and  the  two  pieces  of  the  Dirac  mass  term.  The  mass  Lagrangian  for  matter 
fermions  is  thus 


(4.7) 

The  mass  matrices  are  given  below  in  terms  of  the  parameters  in  the  su- 
perpotential and  the  veto's  of  the  various  doublets.  It  is  clear  that  the  first  four 
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terms  in  (4.7)  are  Dirac  masses  of  the  form  +  ^-c-  m  contrast,  the  last 

term  is  a  "right-handed"  Majorana  mass  of  the  form  ipR^L- 
The  quark  mass  matrices  are 


0        a  (H2)  0 
M  :   =  4   I  f6(L2)        0  {4.8) 

0  0  b(H2) 


\/2 


/    0       d(K2)        0  \ 
M2  =  -4p     d(K2)    e(K'2)    f(K»)  (4.9) 

V  0    /(*?>  »(^)/ 

The  parameter  p  represents  an  S£/(3)  renormalization  group  enhancement 
factor  between  the  GUT  scale  and  low  accelerator  energies.  The  parameter  p 
is  determined  by  any  one  well-known  ratio  of  eigenvalues  of  M_i  and  Mi  at 

3 

low  energies  (see  below),  neglecting  differences  in  the  SU(2)  and  U(l)  renor- 
malization of  the  fermion  masses. 
The  lepton  mass  matrices  are 


3c 


(L2)       0  (4.10) 


0       d(K2)  0 
d(K2)    e(K'2)    f(K»)  |  (4.11) 
^    0       i{K«)  g(K2) 

(    0      ixMo      0  \ 
Mr  =     iiM0      0         0  (4.12) 
V    0         0  j2M0J 

Now  the  two  doublets  from  H  and  L  (coupling  to  the  5  x  10  sectors)  will  mix 
with  doublets  from  other  heavy  fields,  leaving  light  only  the  linear  combination 
cos  a  H2  +  sinaL2-  Similarly,  the  only  doublet  from  higgses  coupling  to  the 
10  x  10  sector  which  remains  light  is 


axK  +  a2K'2  +  a3K'2'  (4.13) 
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Inverting  the  unitary  mixing  matrices,  we  can  express  the  ueu's  in  the 
mass  matrices  (4.8)-(4.12)  in  terms  of  the  vev's  afp  and  oj<t  °f  the  light  higgs 
doublets: 

(#2)  —  °l  =  cosaofj* 
(L2)  =  02  =  —  sin  a  o pj> 

(K2)  =  a3  =  a{  aTT  (4.14) 
(K J)  =  a4  =  °TT 
{K'i)  =  a5  =  a*  arr 

with 

El".|2  =  "2  =  =  (188  GeV)2  (4.15) 

i=l  9 
We  have  enough  freedom  to  redefine  the  righjt-  and  left-handed  fields  so  that 

we  can  make  the  e*j  real  by  rotating  their  phases  away. 

Proceed  now  to  find  the  eigenvalues  of  the  above  mass  matrices  and  the 

mass  eigenstates      as  a  linear  combination  for  the  current  eigenstates  which 

we  have  dealt  with  so  far, 

V>  =  Vxjjm  (4.16) 

with  V  unitary. 

For  charged  fermions,  all  mass  terms  are  Dirac,  of  the  form  (+h.c.) 

tRMxl>L  =  4>RVV-lMVV-lj>L  =  VwAVw,  (4.17) 

where  A  =  V~1MV  is  diagonal.  These  mass  eigenvalues  can  be  made  real 
and  positive  by  further  rotating  the  charge-conjugate  ("right-handed")  fields 
by  some  phases.  Indeed,  right-handed  and  left-handed  fields  mix  differently  in 
general. 
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For  charge-(-j)  fermions  gettting  their  mass  from  Of?, 

(cos  $i     sin  6 1    0  \ 
-sin^    cos0£    0  (4.18) 
0  0  1/ 


and 


so  that 


(cos  0C      sin  0C    0  \ 
-sin0c    cos0c    0  (4.19) 
0  0  1/ 


and 


with 


=  diag(—me,mn,mT) 

..    I    1  (qgl)2    3  1 

A_i   =  diag(md,  -m3,mb) 

3 


cos  20c  =   -,  i.e.,  de- 
nts +  m<i 


and 


mM-me    .  /m7 

cos  20*  =  — -  ,  i.e.,  0£  ~  .  /  

mM  +  me  Y  mM 

The  Georgi-Jarlsk0g  relationships19  are 

m^x  9  ms 
rrifi  ms 


(4.20) 


(4.21) 


(4.22) 


(4.23) 


(4.24) 


=  3—  (4.25) 

and  we  can  identify  the  Yukawa  couplings  a,  b,  c,  as  well  as  the  renormalization 
factor  p: 


=  —y/2m 


-  ^-mT  (4.26) 


2 


p  =  ^  =  2.81  (4.27) 
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To  make  masses  positive,  let 


Va  =  Vt  ,    V_kL  =  V_!  (4.28) 

3  3 

and  redefine  the  charge-conjugate  fields  such  that 

Va  =  Vt^    1    J  ,    V_,R  =  -1    J  (4.29) 

For  the  charge- §  quarks,  which  get  their  masses  from  ayj-,  first  disregard 

entries  in  M2  proportional  to  d  (<<  e,f,g),  which  means  disregard  the  up- 

3 

quark  mass.  There  is  then  enough  freedom  to  absorb  all  phases  in  M2  by 

3 

redefining  the  c-  and  i-quark  fields  (and  their  charge-conjugates).  Assume, 
equivalently,  that  ea4,  do^,  go$  are  real.  Then 


with 


so  that 


10  0 

=   (  0     cost;     sinr;      +  0(  — )  (4.30) 
0    —sinr/    cosrj/  mc 


r/  — —        <<  1  (4.31) 

0ct3  -  ea4 


A 2  diag[0,mc,mt] 


3 

-  4p  diag\0,  ea4  -  2r}fa5,ga^  +  2r//a5] 
The  Yukawa  couplings  are  thus 

-mc 

e  =  

4pa4 

r?(m(  -  mc)  <r4 

4pa5  o5 
tut  +  mc 

9  =  - 


(4.32) 


(4.33) 


4/0<73 

Reintroducing  now  the  entry  da$,  complex  with  phase  j£,  we  can  find  CP 

violation  in  the  full  mixing  matrix  V2.  First  examine  the  determinant  of  the 

3 

mass  matrix 

detM?.  =  -elS  mumcmt  =  -{da^)2(gcr3)  (4-34) 
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where  etS  is  the  determinant  of  VgL.  The  full  mixing  matrix  is  then 


V2  * 

3 


f   cosOu        sin0ue*2  —T)s'mOu^el'< 
-sin0u      cos^e^       ~11  +  ^e~t^ 

1 


,4\ 


V 


-rj  sin  0U    -^  +  r/e,-5 


where 


and 


sin0u  ~  0U  ~  ./— 
V  mc 


A  2  =  <iia<7[— mu,  mc,  mt] 

3 


with  mc  and  mt  given  by  (4.32)  and 

mu  =  —lp  

e<74 


so  that 


\do%\  =  —yjmumc 
4p 


Letting 


V?r  =  V2  ,    V2R  =  V2 

3L  5  3H  3 


and  recalling  that  the  hadronic  weak  current  is 


(4.35) 


(4.36) 


(4.37) 


(4.38) 


(4.39) 


(4.40) 


(4.41) 


one  can 


check  that  the  Kobayashi-Maskawa  matrix,62  KM,  ends  up  being 


3R  3 

COS  0C 


y.\jo  vc  sin  0C  0 

—  sin  0C  cos  rj      cos  0C  cos  77      sin  r] 
sin  0C  sin  r\      —  cos  0C  sin  rj    cos  r/ , 


(4.42) 


up 

(4.34) 


sin  0C   Slfl  77  v. v(j   0111  ry  v-wo 

to  CP-violating  effects,  better  studied  with  the  help  of  the  determinant 
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KM 


u 
c 

t 


V 


(4.43) 


J 


This  prediction  is  to  be  compared  with  the  KM  matrix  obtained  from  an 

analysis  of  experimental  data,  which  yields63 

d  s  b 

0.9733  ±  0.0024    0.231  ±  0.003        <  0.0052 

0.24  ±  0.03        0.97  ±  0.007  0.0435  ±  0.003 
~  0                <  0.076  ~  0.999 

The  best  fit  of  (4.42)  to  (4.43)  provides11 

6C  =  (13.33  ±0.16)° 

t]  =  (2.49  ±0.01)° 

so  that  we  can  establish  (in  disagreement  with  the  usual  Fritzsch  mass  matrix 
result,64,65  and  believing  the  top  mass  to  lie21  around  45  GeV)  that 

mc 


(4.44) 


mt 


(4.45) 


We  can  now  bound  the  Yukawa  couplings  involved,  because  it  follows  from 
(4.15)  that  |at|  <  188  GeV.  We  use  the  values  of  p  and  tj  given  by  (4.27) 
and  (4.31),  the  standard  masses  for  charged  leptons,  and  masses66  of  (5  Mev, 
1.84  GeV,  45t  GeV),  where  t  is  a  dimensionless  parameter  assumed  by  certain 
wishful  propagandists  close  to  one. 


\d 

W.I/ 

\9 


>  5.5  x  10~5 


>  1.3  x  10 


-2 


>  6.5  x  10 


-4 


,-5 


(4.46) 


>  4.5  x  10" 

>  8.7  x  10"4 

>  2.1  x  10"2*2 

To  find  the  eigenvalues  and  eigenmasses  for  neutrinos  (  left-handed  ones 
are  naturally  light  due  to  the  see-saw  mechanism67'68),  we  must  diagonalize 
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a  6  x  6  matrix  which  takes  a  particular  form,  with  the  upper-left  3x3  block 
equal  to  zero.  Following  the  construction  below,69  we  can  reduce  the  problem 
to  the  diagonalization  of  two  3x3  matrices  instead. 

In  the  basis  {ue,!/^,^,!/^^^,^),  the  full  mass  matrix  is 


There  exist 


H  =  (  °,    I)  (4.48) 


and 

with  U\  and  tAj  unitary,  such  that 

VT*>V  =  (%L    A°J  (<•«>) 
with  A/,  and  A#  diagonal,  real  and  positive,  with 

V  =  eiHU  (4.51) 

Indeed,  let 

S  =  -iM*D{M*R)-1  (4.52) 
and  let  U\  and  i/2  be  the  unitary  matrices  such  that 


mi  =  -MDM~lMl 


and 


A*  =  U2Tm2U2 


(4.53) 


1  +  1    ^  (4-54) 

™2  =       +  -(Af^)~  MpAfp  +  -MlMiiMi)-1 

One  can  then  predict  neutrino  masses  and  oscillation  lengths  from  the 

parameters  in  the  mass  matrices,  which  are  more  or  less  known  for  Mj)  and 
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open  to  speculation  for  Mr.  The  flavor  of  these  speculations,  rather  than  a 
solid  result,  can  be  obtained  with  a  couple  of  simplifications. 

In  the  model  at  hand,  asume  first  that  d/d=e/e  —  f/f  =  g/g=  £•  This 
assumption  would  follow  from  an  S'O(lO)  symmetry  broken  at  a  higher  scale 
than  the  SU(5)  unification  we  are  considering.  Such  further  unification  must 
be  rejected  in  the  light  of  the  discussion  on  baryogenesis  (see  below).  Let  xt  = 
rr1  =  (jj-Mn)-1,  and  let  a  be  any  of  the  a,-. 

The  eigenmasses  are  then,  in  terms  of  the  quark  masses, 


muJmumc  2    mcJmumc  2    m2  2 

 %  C  ,   Y  6  i  — L£  i  ri  +  ei,  ^2  +  ^2 

2rj  ri  r2 

with 

KO2 


2ri 

and 


(4.55) 


(4.56) 


(2  =  (4.57) 
r2 

The  90%  confidence-level  limits  on  neutrino  masses70'71  improve  with  every 
new  analysis  of  e+e~  runs  and  yet,  given  their  order  of  magnitude 

mVfl  <  250  KeV 

(4.58) 

mVr  <  125  MeV 

they  provide  only  very  weak  constraints.  Tighter  ones  can  be  obtained  from 
Zeldovich's  cosmological  limit72,73 

mVe  +  mU(t  +  m„T  <  100  eV  (4.59) 

implying 

rx  <  104GeV 

(4.60) 

r2  <  108GeV 
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so  that 

{M?J  Ml) 

where  M0  is  measured  in  GeV  (109  <  M0  <  1012).59 

The  mixing  scheme  depends  very  heavily  on  the  relative  magnitudes  of  the 

various  parameters.  Let  us  assume  that      =      =  £-1>  and  neglect  d  (i.e., 

neglect  the  up-quark  mass)  from  all  expressions.  Then 

/l     0  0\ 
UX  «      0     1     I  (4.62) 

v°  -1  V 

where  we  have  taken  (74  ~  as  (=>  /  ~  — e).  The  estimate 

V22  «  1-^2*V4)2  (4.63) 

can  be  identified  with  the  cosine  of  the  effective  mixing  angle  for  muon  neutrino 
disappearance, 

aVfL  —  fcre<74  »  —  <  10"8  (4.64) 
The  other  parameter  relevant  for  neutrino  oscillations  is 


AM2  = 


2  2 
m     -  mvx 


2 

2    „  (4.65) 


-  "Vr  - 


<  3  x  106eV2 

where  we  have  used  mt  a  45  GeV.  Given  that  the  bound  on  j\  is  weaker,  the 
mixing  angle  could  be  brought  up  to  about  10-4.  An  increase  in  the  mixing 
angle,  however,  would  probably  imply  an  increase  in  Am2,  so  the  outlook  is 
rather  grim:  neutrino  oscillations  seem  to  be  at  least  two  orders  of  magnitude 
away  from  currently  imaginable  detectability.74 
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The  Full  Superpotential 
In  order  to  spontaneously  break  SU(5)  at  a  scale  M^-,  and  the  anomalous 
Peccei-Quinn  symmetry  at  a  scale  Mq  (109  <  Mo  <  1012GeV,  from  the  exis- 
tence of  red  giants  and  from  axions  not  closing  the  universe59),  we  add  a  term 
to  the  superpotential11, 

G2  =  i?i(EE;-M^3  +  i72(M2 -Af£)04  +  r/3E3  +  774£'3  (4.66) 

The  SU(5)  content  of  the  various  fields  introduced,  as  well  as  their  charges 
under  the  global  C/(l)'s,  are  shown  in  Table  II. 

In  order  to  give  mass  to  color  triplets  while  leaving  the  isodoublets  light,60'61 
we  introduce  the  term11'12 

G3  =  {XiK  +  \2K' +  \3K")MU  +  faK  +  p2K'  +  foK")MU' 

(4.67) 

+  {liH  +  l2L)EV 

which,  taking  M  »  (£)  =  Mx,  (a)  gives  mass  to  two  combinations  of  the 
quintets  K,  K',  K"  pairing  them  with  the  5  fields  U  and  U',  leaving  one  light 
quintet,  and  (b)  gives  mass  to  one  quintet  (or  rather,  a  doublet  and  a  triplet) 
formed  with  L  and  H. 

Adding  furthermore  a  term11'12 
G4  =  (SxHQi  +  62L®i  +  <53J02)E  +  6AM®i®2 

(4.68) 

+  (piKe'2+  P2K'G'2  +  p3K"Q'2  +  ptJ'O'jX'  +  PsMe'^ 
masses  can  be  given  to  the  color-triplet  components  of  the  two  light  quintets 
formed  with  G3,  pairing  them  with  the  color  triplets  in  J  and  J'.  We  thus  end 
up  with  four  triplets  with  mass  around  the  intermediate  scales  6MX/M  and 
pMx/M  (two  each),  and  four  massless  doublets,  which  will  acquire  a  mass 
~  M\y  from  spontaneous  electroweak  symmetry  breaking.  Notice  that  only 
two  of  the  doublets  (triplets)  couple  to  matter  at  the  tree  level.   From  the 
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discussion  on  the  renormalization  group  equations  below,  we  shall  obtain  an 
expression  for  My.  The  upper  and  lower  bounds  on  the  triplet  masses  translate 
into  the  bounds 

1(T5  <Pi,6i  <  1(T4  (4.69) 

The  anomalous  Peccei-Quinn  symmetry  X  is  not  broken  explicitly;  hence, 
dimension  five  operators  are  suppressed.11,12  The  full  matter  superpotential 

G  =  Gx  +  G2  +  G3  +  G4  (4.70) 

is  anomaly-free. 

The  Z3  discrete  symmetry  D  under  which 

E        e¥  E  (4.71) 

prevents  terms  like  KJ'H  or  H JE  from  appearing  in  G4.  It  effectively  forces  E 
to  be  a  75  instead  of  a  24,  thus  allowing  the  missing  multiplet  mechanism  im- 
plemented here  to  work.  Charges  under  D  are  displayed  in  Table  II.  Familiarity 
is  unfortunately  not  preserved  in  the  full  G,  so  it  is  perhaps  useful  to  think 
of  it  as  a  low-energy  quasisymmetry.  No  familons  are  to  be  expected  in  this 
model,  although  one  could  imagine  more  complicated  versions  of  G2  +  G3  +  G4 
that  would  require  them  at  some  intermediate  scale. 
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TABLE  II.  Field  content  and  symmetries  (vect- 
orial, Peccei-Quinn,  discrete  Z3)  of  G. 


Superfield 

SU(5)  irrep 

V 

X 

D 

H 

5 

2 

-2 

2 

L 

45 

2 

-2 

2 

K 

5 

-2 

-2 

2 

K' 

5 

-2 

-2 

2 

K" 

5 

-2 

-2 

2 

*l 

1 

-10 

-2 

0 

02 

1 

10 

2 

0 

<t>3 

1 

0 

0 

0 

<t>A 

1 

0 

0 

0 

£ 

75 

0 

0 

1 

E' 

75 

0 

0 

2 

U 

5 

2 

2 

1 

u' 

5 

2 

2 

1 

V 

50 

-2 

2 

0 

J 

5 

-2 

2 

2 

J' 

5 

2 

2 

0 

©1 

50 

-2 

2 

0 

e; 

50 

-2 

-2 

1 

e2 

50 

2 

-2 

0 

e'2 

50 

2 

2 

2 

PROTON  DECAY 

We  know75,76'77  that  all  possible  proton-decay78  operators  are  of  one  of  the 
forms  below.  Greek  indices  run  from  one  to  three  and  represent  color,  Latin 
indices  from  the  middle  of  the  alphabet  run  from  one  to  two  and  represent  weak 
isospin  and  Latin  indices  from  the  beginning  of  the  alphabet  run  from  one  to 
three  and  represent  family  or  generation.  Chirality  is  indicated  for  definiteness: 

°lbcd  =  {^aaLufibR)^icR^ndL)€ap1eij 

°lbcd  =  {?iaaR?j0bL){eCcLuldR)tafatij 

3  ~c  c  (5-l) 

°abcd  =  {daaLupbR){eCcLuldR)^1 

°\bcd  =  tiiaaRQjpbLWmcRQn^dLjtapitijtrnn 
Recall  that,  as  usual,  xfriXR  =  Xi^pRi  and  notice  that 

Olbcd   =   °bacd  (5-2) 

The  form  of  these  operators  is  constrained  by  the  requirement  of  gauge 
invariance  under  the  standard  SU(3)  x  SU(2)  x  U(l).  They  arise  dynamically 
from  integrating  away  heavy  particles  (higgs  fields  or  gauge  bosons  of  broken 
symmetries).  The  effective  Lagrangian  for  the  fermionic  components  of  matter 
fields  contains,  in  our  case,  O1  from  (FT)2  terms  and  O2  from  (TT)2,  with  the 
family  structure  dictated  by  familiarity  [J)  invariance.  There  are  no  O3  or  O4 
operators,  though,  because  we  do  not  allow  H  —  K  nor  L  —  K  mixing  terms; 
i.e.,  we  suppress  proton-decay  operators  involving  higgsinos28'29'30  -with  a 
chirality-flipping  mass  insertion  in  their  propagator-  and  gluinos. 
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To  illustrate  the  point,  consider  a  one-family  toy-model  with  a  Yukawa 
superpotential  of  the  form 

aFTH  +  bTTK  (5.3) 

where  F,  T,  H  and  K  transform  as  5,  10,  5  and  5  of  SU(5).  Both  H  and  K 
contain  color  (anti)triplet,  isosinglet  pieces.  These  couple  both  to  lepto-quark 
and  to  diquark  superfields,  causing  proton  decay  in  two  different  channels.  FF 
can  go  into  TT  via  H  exchange,  and  TT  can  go  to  TT  via  K  exchange. 

Because  of  the  helicity  structure  required  by  the  structure  of  chiral  vertices, 
these  processes  correspond  to  "dimension-six"  proton-decay  operators,  in  which 
all  four  external  lines  are  fermionic  and  the  bosonic  component  of  the  triplet 
(from  H  or  K)  propagates  through  the  internal  line. 

Were  we  to  further  allow  the  mixing  term 

fiHK  (5.4) 

then  tripletinos  (but  not  bosonic  triplets)  would  mix,  inducing  the  process 
squark  +  squark  — ►  antiquark  +  lepton  in  the  channel  FT  — ►  TT.  The  en- 
gineering dimension  of  such  a  diagram,  with  two  scalar  and  two  fermionic 
external  legs,  is  five.  The  initial  state  can  be  virtually  produced  from  a  quark- 
quark  pair  (in  the  nucleon)  via  gluino  exchange.  Since  gluinos  couple  strongly, 
the  amplitude  for  the  whole  process  (gluino  exchange  and  higgsino  exchange) 
is  essentially  the  same  as  for  the  dimension- five  higgsino  exchange  subprocess. 
This  amplitude  crucially  differs  from  the  one  for  quark  +  quark  — ►  antiquark  + 
lepton  via  higgs  boson  exchange,  in  that  it  involves  a  triplet  fermion  propaga- 
tor (rather  than  a  triplet  bosonic  one),  so  it  is  suppressed  by  the  mass  (rather 
than  the  mass  squared)  of  the  triplet  superfield. 
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Proton  decay  thus  proceeds  faster  through  these  higgsino-mediated  dia- 
grams (dimension  five)  than  through  higgs-boson  mediated  ones  (dimension 
six)  roughly  by  a  factor  of  (  wH<g!"  )2.  If  the  theory  contains  dimension  five  op- 
erators,  the  mass  of  the  supertriplet  must  be  boosted  above  ~  1014  GeV  (rather 
than  above  ~  1010  GeV)  in  order  to  prevent  fast  proton  decay  rates.60-79-80'81'82 
We  can  suppress  the  mixing  term  nHK  at  the  cost  of  losing  the  usual  baryo- 
genesis  mechanism,  from  the  interference  between  H  — ►  TT  at  tree  level  after 
/x-mixing  of  H  into  K,  and  H  —>  TT  at  one  loop  with  T,  F  and  H  -  K  in  the 
internal  propagators.  Note  that  in  this  one-family  toy  model  this  interference 
is  zero  even  with  the  mixing  term  because  the  product  of  Yukawa  couplings 
is  necessarily  real.83  (The  chiral  structure  of  the  processes  is  what  interests 
us  now.)  Baryogenesis  seems  to  require  an  interaction  among,  and  hence  the 
existence  of,  different  families! 

The  baryogenesis  problem  will  be  solved  shortly,  but  we  can  foresee  that 
since  the  higgs  triplets  will  have  to  be  rather  light,  the  dangerous  mixings 
must  be  suppressed.  The  action  presented  above  does  so  naturally  with  two  of 
the  global  1/(1) 's,  Peccei-Quinn  and  familiarity.  Around  the  electroweak  scale, 
nevertheless,  supersymmetry-violating  effects  will  induce  such  mixing,  safely 
suppressed  by  O(^-J^). 

Proton  Decay  Lagrangian 

The  superpotential  arising  from  (Eq.(4.l))  involving  matter  and  color- 
triplet  higgses  is  given  below,  where  family  is  indicated  explicity  but  color  and 
Dirac  indices  have  been  dropped. 

In  four-component  notation,  rf>x  stands  for 
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For  calculational  purposes,  it  is  much  easier  to  use  two-component  notation 
and  the  physically  relevant  left-  and  right-handed  fields,  substituting  in  all  the 
previous  expressions 

tj)  — ► 

(5.6) 


0C   ►  °%$R 

and  transposing  the  first  entry  in  rpx  after  the  substitution: 


GtripUu  =  +  uc2d\  -  tug^  -  t2iqijV) 


(5.7) 


-  4d/f3  (u^  +  u\t\  +  qutojt')  ~  *9K3  {u\e%  +  u3d3) 

-  4eK'3  {u\t\  +  u2d2)  -  AfK'i  (uc2ec3  +  u%t\  +  q2iq^ij) 

where  £  =  (^J,  q  =  ^ ^ ,  all  the  triplets  have  canonical  kinetic  terms,  and 
the  Yukawa  couplings  a,  b,  c,  d,  e,  /,  g  were  given  in  terms  of  measured  fermion 
masses  and  constrained  doublet  vev's  above. 

Triplets  coupling  to  the  same  matter  in  Eq.  (5.7)  mix, 

(H3  \  _  (  cos  a  sin  a  \  /  5 \  /  cosa  5^  H   \ 
L3  J  ~  y  -  sin  a    cos  a  J  \  X  J        \  -  sin  a  3 ft  H  / 


K, 
K: 


K|  i  =  l  4  A  ij 


2  "2  '2 
t     Alt  J 


«2 

JU) 

XT 

•7 

(5.8) 


Only  S/ry  and  Hjt  remain  light,  with  masses  j-,  Mj^  around  the 
intermediate  scale,  light  by  unification  standards.  The  a,  can  be  chosen  to 
be  all  real,  just  as  we  did  with  the  doublets  when  we  were  dealing  with  the 
fermion  masses.  Although  we  have  used  the  same  notation  (o^)  for  the  entries 
of  the  mixing  matrices  of  doublets  and  triplets,  there  is  no  reason  why  they 
should  be  the  same. 


33 


The  superpotential  (5.7)  thus  becomes  effectively 


(5.9) 


This  is  the  expression  from  which  we  will  compute  the  baryon  asymmetry 
of  the  universe.  Since  5 and  Hyj-  do  not  mix,  we  need  not  consider  higgsinos. 
The  interaction  term  in  the  Lagrangian  involving  the  scalar  components  of  the 
higgs  superfields  5  and  the  fermionic  components  of  matter  in  (5.9)  is  exactly 
the  same  expression  (5.9),  with  the  usual  even-i?  fields  instead  of  superfields. 

At  energies  much  lower  than  the  masses  of  Ef  t  or  &TT  we  can  integrate 
away  Sj^y  and  Ej>x-  Keeping  only  terms  with  only  one  lepton  field,  we  arrive 
at  the  following  low-energy  Lagrangian  for  proton  decay,  where  all  fields  are 
still  current  and  not  mass  eigenstates: 

C  =        (Labcd  <  obcd  >v  -Labcd  <  abed  >e  +  Rabcd[abcd})  (5.10) 

The  L's  and  iE's  are  coefficients  with  indices  in  family  space  (a,  b,  c,  d  € 
{1,2,3}).  The  various  proton  decay  operators  are,  in  the  notation  of  (5.1), 


so  we  are  breaking  up  the  SU(2)  doublets  q  and  I  appearing  in  (5.1),  since 
proton  decay  takes  place  at  energies  below  the  SU{2)  breaking  scale. 


°lbcd  =  (°M  +  \bacd\ 


(5.11) 
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(5.12) 


Explicitly, 

<  abed  >„  =  (dcaucby  {ucdd) 

=  (<EL«*Wl)  (KRddL) 

<  abed  >e  =  {dcaucby  {tcud) 

=  (daLHR)  &RudL) 
[abed]  =  (ucadcb)  {ecud)^ 

=  (KndbL)  {^cLudR) 

or,  in  two-component  notation, 

<  abed  >„  =  -  (ufRdaR^  \ufLddLJ 

<  abed  >e  =  -  (u^daj?)  {e^LudLJ  (5.13) 

[a&cd]  =  -  (dJRuaRj  (e1[LudL) 


TABLE  III.  Coefficients  of  baryon-number-violating 
four-Fermi  operators  in  effective  Lagrangian  before 
electroweak  breaking. 


abed 

2MpTLabcd 

1212 

i  i 
—a  cos  a 

d2a\ 

1222 

ac  sin  a  cos  a 

de  <x\a2 

1223 

dfctictz 

1233 

—ab  cos2  a 

dga\ 

2222 

i  -  i 
—e  sin  a 

e2a2 

2223 

2233 

6c  sin  a  cos  a 

ge  a.\<X2 

2323 

f24 

2333 

9f<*ict3 

3333 

-62  cos2  a 

g2a\ 
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The  coefficients  and  Rabcd  satisfy  the  identities 

Labcd  =  L(ab){cd)   ~  L(cd)(ab) 

(5.14) 

Rabcd  =  R{ab){cd)   ~  R(cd)(ab) 
and  they  are  all  zero  except  the  ones  in  Table  III  and  their  symmetrizations 
according  to  (5.14). 

Rotate  now  all  the  quark  and  lepton  fields  into  the  mass  eigenstates,  de- 
noted by  the  additional  subscript  ra,  and  related  to  the  current  eigenstates  by  a 
unitary  transformation  ^  =  V^VVn-  The  new  proton  decay  operators  have  the 
same  form  as  the  old  ones,  but  with  massive  fields  in  them.  They  are  related 
to  the  old  ones  by  four  mixings,  one  for  each  field.  Summing  over  the  various 
neutrino  flavors,  indistinguishable  in  proton-decay  experiments,  and  neglect- 
ing the  see-saw  coupling  between  ordinary  neutrinos  and  massive  right-handed 
ones,  we  can  express  the  Lagrangian  (5.10)  in  terms  of  mass  eigenstates: 

£  =  £  (Kbd  <abd>  +  Kbcd  <  abcd  >  +Kbcd\abcd})  (5-15) 

The  coefficients 

Kbd  =  £  W^^ae(^)/(^^ 

Kbcd  =  ~    £  Lefgh iV! I R}ea{Vl^)\ [VIl)9c {V- l^h  (5.16) 
e,f,9,h 

Kbcd  =  £  ^m^iji)*-^-^)^^)/^)/ 

e,f,g,h 

can  be  evaluated  using  symbolic  manipulations  on  a  computer,  without  much 
problem.  The  full  formal  expressions  are  long  and  unilluminating. 

Renormalization  of  Proton  Decay  Operators 
One  last  step  must  be  taken  to  obtain  the  effective  proton  decay  La- 
grangian, useful  to  compute  proton  decay  rates:  we  must  renormalize84  all 
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the  (dimension  six)  operators  down  from  unification  to  our  energies.  1  • 
Indeed,  the  expressions  for  the  coefficients  Lv ,  Le,  Re  of  the  various  operators 
are  valid  at  the  unification  scale.  The  renormalization  of  the  Yukawa  cou- 
plings themselves  is  taken  into  account  by  using  low-energy  fermion  masses, 
but  the  full  operators  undergo  further  renormalization  at  low-energies,  by  a 
multiplicative  enhancement  factor  which  can  be  calculated  using  dimensional 
regularization.  We  neglect  U(l)  renormalization,  assume  exact  supersymmetry 
above  M\y,  and  renormalize  with  one-loop  equations  down  to  the  charm  mass. 
Below  M\\r,  we  consider  only  SU(3)  renormalization.  These  are  all  reasonable 
assumptions  in  light  of  the  requisite  accuracy. 

Both  operators  O1  (mediated  by  scalar  #3  and  L3)  and  O2  (mediated  by 
scalar  K$,  K'3,  K'£)  could  mix  under  SU(3)  renormalization  with  their  super- 
symmetric  partners,  the  dimension  five  operators  mediated  by  higgsinos  with 
the  same  quantum  numbers-were  those  operators  allowed  by  the  theory  at  all. 
They  cannot,  however,  mix  with  any  other  operators  under  SU(2)  renormal- 
ization because  under  exact  supersymmetry  the  SU(2)  gaugino  is  Majorana 
massless:  were  it  to  be  exchanged  between  two  external  legs  with  SU(2)  quan- 
tum numbers,  it  would  undergo  a  helicity  flip  to  preserve  the  chiral  structure 
of  the  vertices,  which  it  cannot  do  because  (as  advertised)  it  is  massless.85 

Let  us  consider  the  spectrum11  contributing  to  the  renormalization-group 
equations  between  Mx  (unification  scale)  and  the  low-energy  world.  We  have 
three  light  fermion  families  (the  singlets  are  heavy  but  irrelevant  for  group 
behavior),  with  the  top-quark  mass  around  45GeV,  four  light  isodoublets, 
with  masses  around  M\y,  and  four  light  color  triplets  with  masses  around 
M0  a  1010GeV.  We  assume  all  other  higgs  superfields  have  masses  of  at  least 
0(MX). 
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If  the  coupling  constant  for  the  ith  simple  factor  of  SU(5)  is,  to  one  loop, 
of  the  form 

Q'(M)  =  wdpffl  (5,I7) 

we  can  find  the  coefficients  6^  from  the  following  expressions,  where  nq  (n.£,  nd, 
nt)  is  the  number  of  quark  flavors  (lepton  flavors,  isodoublets,  color  triplets) 
below  the  energy  scale  //, 

For  nonsupersymmetric  QCD,  (SU(3)  below  M\y), 

63  =  11  -  (5.18) 

For  supersymmetric  QCD,  (5(7(3)  above  M\y), 

63  =  3  x  3-  \{2nq  +  nt)  =  3  -  \nt  (5.19) 

where  the  2  comes  from  counting  quarks  and  antiquarks. 
For  supersymmetric  SU(2), 

62  =  3x2-^  +  ^  +  ^)  =  -^nd  (5.20) 

Finally,  for  supersymmetric  U(l), 

1  3 
h  =  --{nq  +  ne)-  frac310nd  =  -6  -  —nd  (5.21) 

As  boundary  conditions  for  the  renormalization-group  equations  we  use 
«7m{Mw)  =  128.5 
a^{Mw)  =  0.12-1  =  8.3(±0.4) 


(5.22) 


Recall  also  that 


a2  1{M\y)  =  aer*(A%)  sin2  6W 


w 


3-1  (5'23) 
al  {M\y)  =  g "cm iMw)  cos2  0W 

3 


(5.25) 
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whereas  the  unification  scale  is  defined  by 

a  =  a^Mx)  =  a2{Mx)  =  a3{Mx)  (5.24) 
The  one-loop  renormalization  equations  are  thus 

.--fdM  =  s<-0*(s£) 

Solving  them,  one  easily  finds11'86 

Mx  =  5.6  x  1016  GeV 

a  =  a{Mx)  =  (19.60)-1  (5.26) 
xw  =  s'm2ew(Mw)  =  0.237 

with  errors  of  about  5%  arising  mostly  from  uncertainites  in  A-jg^,  where  MS 
is  the  modified  minimal  substraction  scheme  in  QCD.  The  rather  high  value 
for  xw  is  in  excellent  agreement  with  the  most  recent  measurements  of  xw  in 
deep-inelastic  neutrino  scattering.87'88 

Now  for  the  actual  renormalization  of  the  qqql  operators.  Since  we  are 
neglecting  U(l)  effects,  all  operators  at  hand  get  renormalized  in  the  same  way. 
Letting  A  be  any  of  the  coefficients  Lu,  Le,  Re  in  the  effective  Lagrangian,  the 
enhancement  factors  relate  A  at  one  scale  to  A  (and  all  other  B  with  the  same 
quantum  numbers)  at  some  other  scale.  Specifically, 

A{mp)  ~  A{mc)  ~  E2E3ngE333A(Mx)  (5.27) 
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where  we  have  neglected  E\  and  the  enhancement  factors  are 


x 


<xzimb). 


t57*14' 


(5.28) 


(4) 


The  7's  are  the  relevant  eigenvalues77  of  the  anomalous  dimension  matrices, 
which  happen  to  be  the  same  as  those  recently  considered  for  the  case  of 
dimension  five  operators:85 


72  =2 

73na  =  2 
4 

88 

73      =  o 


(5.29) 


The  /^-function  coefficients  can  be  read  of  from  (5.18),  with  =  4.  The 
superindex  in  E$ns  denotes  nq,  and  fit  in  ^3S5- 

With  a  top-quark  mass  of  45GeV  [and  a3(M0)  m  17_1,  M0  a  1010  5GeV], 
the  enhancement  factor  is 


E  =  £2  x  £3™  x  E3ss  =  1.39  x  1.26  x  1.65  =  2.89  ±  0.2  (5.30) 


Proton  Decay  Rates 
Armed  with  this  low-energy  potential,  we  can  calculate  the  proton  (and 
neutron)  decay  amplituyes  into  various  channels.  However,  a  nucleon  at  rest  is 
hardly  a  system  made  up  of  three  bound  pointlike  quarks,  and  although  much 
effort  has  been  made,  using  various  techniques,78  to  produce  reliable  ampli- 
tudes in  terms  of  the  qqqi  operators,  large  uncertainties  remain.  We  follow 
here  the  most  straightforward  approach.89  For  the  phase  space,  we  assume  the 
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quarks  to  be  non-relativistic  and  the  antilepton  extremely  relativistic.90  All  nu- 
clear effects  are  disregarded,  and  we  concentrate  on  two-body  p  ->  ml  decays; 
that  is,  we  ignore  the  pion  pole  contributions  and  all  three-body  decays.  Pro- 
ton and  neutron  wave  functions  are  taken  to  be  51/(6)  symmetric,91'92  where 
SU{6)  D  SU{2)spin  x  St/(3)flavor  is  the  light-cone  spin-flavor  symmetry. 

Following  Ref.  89,  the  lifetimes  of  each  mode  can  be  conveniently  written 
in  the  generic  form11 

(5.31) 


r  =  y 


m 


LlOiOJ 

where  m  is  the  mass  (in  GeV)  of  the  color  triplet  mediating  the  decay,  MpT 
if  the  outgoing  antilepton  is  right-handed,  and  M?t  »*  is  left-handed.  The  y 
coefficients  (in  years)  for  the  main  decay  modes  are  shown  in  Table  IV.  Their 
explicit  dependence  on  the  triplet  mixing  parameters  is  also  shown,  and  since 
we  are  looking  for  an  limit  on  MFt  and  Mjr,  the  upper  bound  (188  GeV)  for 
all  doublet  veu's  has  been  used,  i.e.,  the  bounds  in  (4.46)  have  been  saturated. 

At  90%  confidence  level,  the  latest  limits  on  the  proton  lifetime  into  strange 
modes  are6 

r  (p  ^  fJ.+K°)  >  4  x  1031years 
\  >  (5.32) 

t  (p  -»  uK+)  >  5  x  1031  years 
These  limits  can  be  translated  into  lower  bounds  on  the  masses  of  the  color 

triplets  mediating  proton  decay: 

mFT  ^  0.75  x  Y/I15C2-2.8S2  -0.0555CI  x  1010  GeV 

mFT  >  0.53  x  ^/|3.2c2  -  12.2s2  -  20sc|  x  1010  GeV 

mTT  >  0.32  x  1010  GeV 

x      -  319q2  +  319a2,  +  0.6a|  -  2.7aia2  -  0.17aia3  -  27.8a2a3| 

(5.33) 
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where  c,  s  stand  for  cosa,  sina,  and  cq,  a2,  and  a3  stand,  respectively,  for 
cos  a',  sin  a'  cos  a",  and  sin  a'  sin  a".  For  any  values  of  the  mixing  angles,  the 
above  expressions  can  only  vary  so  much,  so  that  a  perhaps  more  useful  result 
is 

m  fp  = 


(4-10) 


4  x  1031years. 


i 

x  1010GeV 


mj"T 


(5-18)  TvK 


i 

x  10luGeV 


(5.34) 


4  10 


5  x  1031years 

where  the  lifetimes  are  normalized  by  the  current  limits,  for  years  of  use. 
It  follows  that  masses  of  the  color  triplets  are  bound  below  by 

mFT,mTT  >  2xl010GeV  (5.35) 

This  lower  bound  on  the  triplet  masses  translates  into  a  lower  bound  for 
the  mass  of  the  infiaton  {A2/M)  from  which  triplets  decay,7'8'9-10'11'12'13  so 
we  can  extract  a  bound  on  A  purely  from  underground  detector  results: 

A  >  1CT4-25  (5.36) 

Remarkably,  proton  decay  experiments  set  a  lower  limit  on  A  which  is 
essentially  the  same  as  the  lower  limit  (3.9).  The  upper  limit  (3.5)  on  A 
arising  from  cosmological  considerations  on  gravitino  abundance  at  the  time  of 
nucleosynthesis  is  very  close,  so  A  is  (as  repeatedly  advertised  all  along)  very 
tightly  constrained. 

It  is  perhaps  amusing  to  extract  more  mileage  from  the  relationship  be- 
tween the  reheat  temperature  and  the  gravitino  mass,  Eq.(3.4).  The  bound 
(5.35)  translates  into  the  bound  (5.36)  because  inflatons  must  decay  into  triplets 
in  order  to  produce  the  baryon  asymmetric  universe  we  live  in.  The  two  in- 
equalities (3.4)  and  (5.36)  result  in 

H  <  1.54  xlO-7  =  3.7xlOnGeV  (5.37) 
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Proton  decay  experiments  have  thus  obtained  an  upper  bound  on  the  gravitino 
mass!  In  numbers  the  result  is  not  so  impressive: 

m2  <  57TeV  (5.38) 

2 

Along  the  same  lines,  the  lack  of  experimental  confirmation  for  supersym- 
metry  in  accelerator  experiments  allows  us  to  expect 

m3>100GeV  (5.39) 

2 

from  which  it  follows  that 

A  <  1(T3-3  (5.40) 

hence 

ro3  <  1.9  x  1012GeV  =>  t^k  <  1.3  x  1040years  (5.41) 

On  the  other  hand,  assuming  supersymmetry  is  not  found  at  the  SSC, 
then  the  gravitino  mass  will  be  narrowly  constrained  between  lower  bounds 
from  the  SSC  and  upper  bounds  from  upgraded  proton  decay  experiments.  If 
the  largest  conceivable  Earth-based  detectors6  are  built  and  work,  the  proton 
lifetime  limit  will  improve  to 

TnK  >  5  x  1034years  (5.42) 

Then 

A  >  10"3-9  (5.43) 

and 

ms  <  1  TeV  (5.44) 

2 

which  could  be  a  useful  constraint. 
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TABLE  IV.  Decay  rates  for  two-body  final  states 


Mode  T  :  Branching  Rates  (yr  *) 


4-° 

1.02  x  io27|£  x  Lenn\2 

4*° 

1.02  x  1027\E  x  i2iU1|2 

4*° 

1.39  x  1027|£  x  Le2Ul\2 

4*° 

1.39  x  1027|£  x  L\2n\2 

4^° 

1.00  x  1027|£  x  L\ni\2 

4-° 

i.oo  x  io27|£  x  flfmr 

1.34  x  1027|E  x  L|121|2 

1.34  x  1027|£7  x  L\22l\2 

2.05  x  1027|£7  x  L\n\2 

uRK+ 

1.39  x  1027|£  x  Lvn2\2 

VRK*+ 

0.012  x  1027|£  x  (3L^12  -  2L^n)|2 

e+rj 

0.20  x  r[e+7r°] 

e+p° 

0.057  x  r[c+7r°] 

e+u 

0.49  x  r[e+7T°] 

e+K*° 

0  0088  x  Tle+.K'0! 

0.20  X  r[/X+7T°] 

0.045  x  I>+7T°] 

0.38  x  I>+7r°] 

I7p+ 

0.20  x  T\u*+\ 

BARYOGENESIS 
The  observed  baryon  asymmetry2,93  of  the  universe, 

"B-b/""!  ~  10"10  l6'1) 

is  a  quantification  of  two  surprising  facts.  One,  that  no  significant  amounts 
of  antimatter  exist.  Two,  that  in  fact  there  is  much  more  matter  per  unit 
entropy  (per  photon)  than  one  would  expect  from  an  equilibrium  evolution  of 
the  primeval  "singularity."  The  generation  of  this  baryon  asymmetry,  a  process 
known  as  baryogenesis,  requires  three  fundamental  conditions.94  First,  baryon 
number  must  obvioulsy  be  violated.  Next,  thermal  equilibrium  must  also  be 
violated.  Finally,  C  and  CP  must  also  not  be  conserved. 

Baryogenesis  may  proceed  after  the  inflationary  phase  of  the  evolution  fo 
the  universe,  by  allowing  the  matter  sector  of  the  superpotential  G  to  con- 
tain higgs  triplets  which  produce  a  non-zero  baryon-number  asymmetry  SB 
per  triplet-antitriplet  decay,  and  such  that  the  inflaton  can  decay  into  them 
preferentially  (i.e.,  m$  <  m^). 

Since  the  inflaton  is  never  in  equilibrium  with  other  fields,  we  do  not  have 
to  worry  about  producing  enough  triplets  via  a  Boltzmann  distribution.24  In- 
deed, the  inflaton  oscillates  into  the  heaviest  fields  around,  and  we  can  simply 
estimate7'8'11,12  that  all  of  the  inflaton's  energy  is  released  into  heavy  triplets 
which  quickly  decay  into  radiation,  reheating  the  bubble  to  a  temperature  Tr: 
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*  6B^  *  SB^ 

n~f 

~  6B-^=-  *  6B^  (6.2) 

~  SB— 
M 

In  this  expression  S B  is  the  baryon  number  asymmetry  produced  per  decay 
of  an  inflaton  and  an  anti-infiaton.  Given  that  the  inflaton's  branching  ratio 
into  color  triplets  is  practically  one,  then  6B  is  essentially  the  baryon-number 
asymmetry  produced  per  triplet-antitriplet  decay,  which  we  calculate  from  the 
matter  sector  of  the  GUT  superpotential.  Clearly,  G  must  contain  triplets  of 
mass  a  little  lower  than  the  inflaton's,  with  a  rather  high  decay  asymmetry 

SB  ~  10~6  (6.3) 

At  tree  level,  all  |amplitudes|2  are  real,  so  CP  violation  can  be  produced 
perturbatively  only  at  the  one-loop  level.  We  need  one-loop  decays  which  can 
interfere  with  tree-level  decays  to  produce  a  net  baryon  asymmetry/0^0'^''5'0^ 

Indeed,  if  a  generic  particle  X  (colored  higgs,  say)  can  decay  into  two 

channels  with  different  baryon  number,  for  instance 

2 

X  — ►  qq        (-Bfinai  =  --), 

!  3  (6-4) 
X  — y  ql       (-Bfinal  =  -), 

where  q  is  a  quark  and  I  a  lepton  (superfields) ,  then  the  baryon  asymmetry 
SB  produced  per  decay  of  a  pair  of  X  and  X  is 


=  jT(X  -»  ql)  -  |r(X  ->  qq)        -^(X  -»  ql)  +  jfjX  -»  qq) 
T{X-+qe)  +  r{X-+qq)     +      r{X-+qi)  +  r{X-+qq) 
T(X->qq)-T(X^qq)  (65) 

r(X  -+  anything) 
T{X  -*  ql)  -  T(X  ->  ql) 
T{X  ->  anything) 
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where  we  have  used  CPT  invariance. 
In  general, 

sb  =  ^^Bar.-r.)  (e.e) 

where  the  sum  extends  over  all  decay  channels  t  of  X,  each  with  a  net  baryon 
number  B{. 
Let 

T(X -><?£)  =  \gTfT  +  9LfL\2  (6-7) 

where  the  amplitude  is  the  sum  of  a  tree-level  graph  and  a  one-loop  graph 
(disregarding  higher  orders  in  h),  with  g  the  (product  of)  coupling  constant (s) 
and  /  the  dynamical  factor,  which  depends  on  masses  and  external  momenta 
and  is  the  same,  by  CPT  invariance,  for  X  as  for  X  decays.  Then 

T(X-+qt)  =  \g*TfT  +  ffl/ll2  (6-8) 

so  that  (6.5)  becomes  simply 

gB  =  4lm(/^/£)Im(g^L) 
T(X  -*  anything) 

or,  in  the  general  case  with  more  than  two  decay  channels, 


XD  4EBtIm(/^./£t)Im(4tg£t) 

OD    =   ~2  

E  \9Tihi  +9LifLi\ 


(6.10) 


„  ^Bjlm(fyLi)lm(g*TigLi) 
T,\9Tifc\2 

The  expression  lm(fj,fi1)  can  be  calculated  using  the  Landau-Cutkosky 
rules  for  the  interference  diagram  f^lL-  Since  at  tree  level  all  the  phases 
disappear,  we  may  choose  couplings  so  that  fj.  =  fx-  Then  ha(f^fi)  = 
/rlm/i-  To  evaluate  this,  all  one  has  to  do  is  modify  the  internal  propagators 
in  fi  separating  initial  from  final  states  by  putting  them  on  shell. 
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All  said  and  done,  only  the  innermost  line  in  the  diagrams,  where  a  higgs 
triplet  is  usually  propagated,  is  allowed  to  be  off-shell.  We  could  add  gauge 
interactions  to  these  diagrams  but  (contrary  to  low-energy  processes,  such  as 
proton  decay)  at  the  energies  relevant  for  our  dicussion,  all  gauge  couplings  are 
weak,  so  that  renormalization  of  these  operators  can  be  neglected.  All  vertices 
are  chiral,  and  even  cubic  if  we  are  willing  to  introduce  spurious  fields  to  take 
into  account  mass  terms. 

The  only  way  arrows  (chirality)  can  arrange  themselves  in  the  interference 
diagrams  is  with  mass  insertions  (chirality  flips)  in  either  the  external  line  only, 
or  else  in  two  of  the  internal  non-higgs  lines.  The  origin  of  these  constraints  is, 
in  supergraph  language,  the  fact  that  interactions  in  the  superpotential  involve 
only  chiral  fields,  so  all  three  arrows  point  either  into  or  out  of  any  given  vertex. 

All  the  baryon  asymmetry  is  generated  from  the  higgs  fields  coupling  to 
the  10  •  10  sector,  for  otherwise  the  product  of  couplings  g*TgL  is  real.  More 
particularly,  SB  #  0  arises  from  the  interplay  between  the  only  two  higgs 
superfields  coupling  to  the  same  sector,  both  with  non-diagonal  (in  family 
space)  couplings.  Remarkably,  both  higgs  bosons  and  higgsinos  can  produce  a 
baryon  asymmetry  upon  decay. 

Taking  into  account  that  K'J  =  ol\EtTi  Ks  -  ajarr,  one  can  evaluate  the 
interference  and  obtain  the  following  asymmetry  per  decay:11'12 

SB  _   i  |«ii2|a3l2imKg/r/^Ty2) 

327r2|/1|2  +  |/2|2  +  |di|2  +  |eT  +  |y|2 
Given  that  the  triplet  must  decay  into  neutrinos  (right-handed  ones),  with 
masses  given  by  (4.55),  it  is  clear  that 


M0 


(6.11) 
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so  we  can  write,  assuming  a  hierarchy  of  Yukawa  couplings  and  all  phases  to 
be  of  order  unity, 

„  <  i  MWillAllAl  (6.13) 

32tt2  \g\ 

To  estimate  this,  take  aj  ~  «s  ~  1,  disregard  all  doublet  mixings,  and 
assume  proportionality  between  tilded  and  untilded  Yukawa  couplings  [via 
50(10)].  The  value  thus  obtained, 

SB  <  io-3  V^m'  ^  10-io  (6  14) 

mtaL 

is  much  too  low,  so  the  presumed  5O(10)-inpired  proportionality  between 
tilded  and  untilded  Yukawa  couplings  is  to  be  rejected.  Indeed  taking  all 
Yukawa  couplings  in  (6.13)  to  be  of  order  10-1,  6B  is  about  10~5,  an  order  of 
magnitude  above  the  required  value. 

It  is  thus  rather  easy  to  produce  the  baryon  asymmetry  of  the  universe  from 
triplet  decays  into  right-handed  neutrinos  and  ordinary  quarks.  The  parame- 
ters to  play  with  are  those  in  the  neutrino  Dirac  mass  matrix,  while  making 
the  right-handed  Majorana  masses  as  large  as  allowed,  i.e.,  smaller  than  the 
triplet  mass  Mj<t-  The  rather  large  Yukawa  couplings  between  SU(5)  singlets 
and  quarks  could  be  ruled  out  if  neutrino  oscillations  were  to  be  observed,  as 
pointed  out  above. 

Baryon  number  violation  has  been  searched  for  with  greatest  care  in  un- 
derground detectors,  to  no  avail.  The  rather  strong  limits  on  proton  decay 
rates  provide  us  with  a  lower  limit  on  the  mass  of  the  baryon-number-violating 
triplets.  The  upper  limit  on  this  mass  (assuming  Yukawa  couplings  to  be  per- 
turbative,  i.e.,  smaller  than  one)  sets  an  upper  limit  on  the  lifetime  of  the  pro- 
ton, limit  which  unfortunately  falls  beyond  the  neutrino-induced  background 
on  Earth. 


CONCLUSIONS 

We  have  achieved  to  construct  a  model  in  which  the  two  bounds  on  the 
mass  of  the  color  triplets  (lower  from  proton-decay  experiments,  upper  from 
the  requirement  of  baryogenesis)  are  met.  Good  relationships  between  fermion 
masses  and  a  reasonable  Kobayashi-Maskawa  matrix  (including  C P  violation) 
have  been  incorporated  into  an  anomaly-free  grand-unified  model. 

We  have  also  established  a  hierarchy  of  masses 

mrightneutrino   <  mtriplet   <  minflaton  <  mO'Raifearton 

at  the  intermediate  scale,  ~  1010  GeV  This  intermediate  scale  arises  from  a 
plurality  of  physical  considerations  associated  with  supersymmetry  breaking, 
inflation,  baryogenesis,  and  the  breaking  of  the  Peccei-Quinn  symmetry. 
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